Prepared for submission to JCAP 



CMB polarization anisotropies from 
cosmological reionization: extension 
to B-modes 



T. Trombetti, a6 C. Burigana ac 

a INAF-IASF Bologna, Via P. Gobetti 101, 1-40129, Bologna, Italy 
fc Dipartimento di Fisica, Universita La Sapienza, P.le A. Moro 2, 1-00185 Roma, Italy 
c Dipartimento di Fisica, Universita degli Studi di Ferrara, Via G. Saragat 1, 1-44100 Ferrara, 
Italy 

E-mail: trombetti@iasfbo.inaf.it, burigana@iasfbo.inaf.it 

Abstract. The accurate understanding of the ionization history of the Universe plays a fun- 
damental role in modern cosmology. It includes a phase of cosmological reionization after 
the standard recombination epoch, possibly associated to the early stages of structure and 
star formation. While the simple "r-parametrization" of the reionization process and, in 
particular, of its imprints on the cosmic microwave background (CMB) anisotropy likely rep- 
resents a sufficiently accurate modelling for the interpretation of current CMB data, a great 
attention has been recently posed on the accurate computation of the reionization signa- 
tures in the CMB for a large variety of astrophysical scenarios and physical processes. This 
work is aimed at a careful characterization of the imprints introduced in the polarization 
anisotropy, with particular attention to the B-modes. We have implemented a modified ver- 
sion of CAMB, the Cosmological Boltzmann code for computing the angular power spectrum 
(APS) of the anisotropies of the CMB, to introduce the hydrogen and helium ionization 
fractions predicted in astrophysical and phenomenological reionization histories, beyond the 
simple r-parametrization. We compared the results obtained for these models for all the 
non- vanishing (in the assumed scenarios) modes of the CMB APS. The amplitude and shape 
of the B-mode APS depends, in particular, on the tensor-to-scalar ratio, r, and on the reion- 
ization history, thus an accurate modeling of the reionization process will have implications 
for the precise determination of r or to set more precise constraints on it through the joint 
analysis of E and B-mode polarization data available in the next future and from a mission 
of next generation. Considering also the limitation from potential residuals of astrophysi- 
cal foregrounds, we discussed the capability of next data to disentangle between different 
reionization scenarios in a wide range of tensor-to-scalar ratios. 
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1 Introduction 

The accurate understanding of the ionization history of the Universe plays a fundamental 
role in modern cosmology. The classical theory of hydrogen recombination for pure baryonic 
cosmological models [1, 2] has been subsequently extended to non-baryonic dark matter 
models [3-5] and recently accurately updated to include also helium recombination in the 
current cosmological scenario (see e.g. [6] and references therein). Various models of the 
subsequent Universe ionization history have been considered to take into account additional 
sources of photon and energy production, possibly associated to the early stages of structure 
and star formation, able to significantly increase the free electron fraction, x e , above the 
residual fraction (~ 10 -3 ) after the standard recombination epoch at z rcc ~ 10 3 . These 
photon and energy production processes associated to the cosmological reionization phase may 
leave imprints in the Cosmic Microwave Background (CMB) providing a crucial "integrated" 
information on the so-called dark and dawn ages, i.e. the epochs before or at the beginning the 
formation of first cosmological structures. For this reason, among the extraordinary results 
achieved by the Wilkinson Microwave Anisotropy Probe (WMAP) mission 1 , the contribution 
to the understanding of the cosmological reionization process has received a great attention. 

This work is aimed at a careful characterization of the imprints introduced in the po- 
larization anisotropy, with particular attention to the B-modes. 

In Sect. 2 we briefly summarize the current observational information coming from 
WMAP on the cosmological reionization and describe its main imprints in the CMB. In 
Sects. 3 and 4 we discuss the main properties of the two classes of reionization models 
considered in this work, astrophysical and phenomenological, respectively. Sect. 5 concerns 
the numerical implementation we carried out to include the considered reionization scenarios 
in our Boltzmann code modified version. The experimental sensitivity of on-going and future 
CMB anisotropy space missions and the limitation coming from astrophysical foregrounds 
are discussed in Sect. 6, while our main results are presented in Sect. 7. Two appendices 
report on some technical details of this work. Finally, in Sect. 8 we draw our conclusions. 

1 http: / /lambda.gsfc. nasa.gov / product /map / current / 
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2 Cosmological reionization 



To first approximation, the beginning of the reionization process is identified by the 
Thomson optical depth, r. The values of r compatible with WMAP 3yr data, possibly com- 
plemented with external data, are typically in the range ~ 0.06 — 0.12 (corresponding to a 
reionization redshift in the range ~ 8.5 — 13.5 for a sudden reionization history), the exact 
interval depending on the considered cosmological model and combination of data sets [7]. 
Subsequent WMAP data releases improved the measure of r, achieving a 68% uncertainty 
of ~ ±0.015 [8-11]. Under various hypotheses (simple ACDM model with six parameters, 
inclusion of curvature and dark energy, of different kinds of isocurvature modes, of neutrino 
properties, of primordial helium mass fraction, or of a reionization width) the best fit of r 
lies in the range ~ 0.086 — 0.089, while allowing for the presence of primordial tensor per- 
turbations or (and) of a running in the power spectrum of primordial perturbations the best 
fit of r lies in the range to ~ 0.091 — 0.092 (0.096). While this simple "r-parametrization" 
of the reionization process and, in particular, of its imprints on the CMB anisotropy likely 
represents a sufficiently accurate modelling for the interpretation of current CMB data, a 
great attention has been recently posed on the accurate computation of the reionization sig- 
natures in the CMB for a large variety of astrophysical scenarios and physical processes (see 
e.g. [12-20]) also in the view of WMAP accumulating data and of forthcoming and future 
experiments beyond WMAP (see [21] for a review). In [22] a detailed study of the impact 
of reionization, and the associated radiative feedback, on galaxy formation and of the cor- 
responding detectable signatures has been presented, focussing on a detailed comparison of 
two well defined alternative prescriptions (suppression and filtering) for the radiative feedback 
mechanism suppressing star formation in small-mass halos, showing that they are consistent 
with a wide set of existing observational data but predict different 21 cm background signals 
accessible to future observations. The corresponding signatures detectable in the CMB have 
been then computed in [23] . 

Different scenarios have been investigated in [24] assuming that structure formation 
and/or extra sources of energy injection in the cosmic plasma can induce a double reioniza- 
tion epoch of the Universe at low (late processes) or high (early processes) redshifts, providing 
suitable analytical representations, called hereafter as " phenomenological reionization his- 
tories". In the late models, hydrogen was typically considered firstly ionized at a higher 
redshift (z ~ 15, mimicking a possible effect by Pop III stars) and then at lower redshifts 
(z ~ 6, mimicking the effect by stars in galaxies), while in the early reionization framework 
the authors hypothesized a peak like reionization induced by energy injection in the cosmic 
plasma at z ^ some x 10 2 . 

In this work we present a detailed analysis of these relatively wide sets of astrophysical 
and phenomenological models and of the signatures they induce in the CMB temperature and 
polarization anisotropies, but the methods described here can be used as guidelines for the 
implementation of any other reionization scenario (see e.g. [25]). In particular, we present 
original computations for the polarization B-mode angular power spectrum (APS) and com- 
pare them with the sensitivity of on-going (Planck 2 , see [26]) and future (assuming COrE 3 , 
see [27], as a reference) CMB space missions. 



2 www.rssd.esa.int/Planck 
3 http://www. core-mission. org/ 
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2.1 Signatures in the CMB 

The cosmological reionization leaves imprints on the CMB depending on the (coupled) 
ionization and thermal history. They can be divided in three categories 4 : i) generation 
of CMB Comptonization and free-free spectral distortions associated to the Inter Galactic 
Medium (IGM) electron temperature increase during the reionization epoch [30-34], ii) sup- 
pression of CMB temperature anisotropics at large multipoles, £, due to photon diffusion, 
and Hi) increasing of the power of CMB polarization and temperature-polarization cross- 
correlation anisotropy at various multipole ranges, mainly depending on the reionization 
epoch, because of the delay of the effective last scattering surface. The imprints on CMB 
anisotropies are mainly dependent on the ionization history while CMB spectral distortions 
strongly depend also on the thermal history. The reionization process mainly influences the 
polarization E and B modes, and the cross-correlation temperature-polarization mode be- 
cause of the linear polarization induced by the Thomson scattering. The effect is typically 
particular prominent at low multipoles £, showing as a bump in the power spectra, otherwise 
missing. 

Through this note we assume a flat ACDM model compatible with WMAP, described 
by matter and cosmological constant (or dark energy) density parameters Q m = 0.24 and 
= 0.76, reduced Hubble constant h = ffo/(100km/s/Mpc) = 0.73, baryon density Qbh 2 = 
0.022, density contrast as = 0.74, and adiabatic scalar perturbations (without running) with 
spectral index n s = 0.95. We adopt a CMB background temperature of 2.725K [35]. 

3 Astrophysical reionization models 

The analysis of Lya absorption in the spectra of the 19 highest redshift Sloan Digital Sky 
Survey (SDSS) Quasar (QSO) shows a strong evolution of the Gunn-Peterson Lya opacity 
at z ~ 6 [36, 37]. The downward revision of the electron scattering optical depth to r = 
0.09 ± 0.03 in the release of the 3yr WMAP data, confirmed by subsequent releases, is 
consistent with "minimal reionization models" which do not require the presence of very 
massive (M > IOOMq) Pop III stars [38]. The above models can be then used to explore the 
effects of reionization on galaxy formation, referred to as "radiative feedback" . 

A semi-analytic model to jointly study cosmic reionization and the thermal history of 
the IGM has been developed in [39]. According to Schneider and collaborators, the semi- 
analytical model developed by Choudhury &; Ferrara, complemented by the additional physics 
introduced in [40], involves: i) the treatment of IGM inhomogeneities by adopting the pro- 
cedure of [41]; ii) the modelling of the IGM treated as a multiphase medium, following the 
thermal and ionization histories of neutral, HII, and Helll regions simultaneously in the pres- 
ence of ionizing photon sources represented by Pop III stars with a standard Salpeter IMF 
extending in the range 1 — 100 M [42], Pop II stars with Z = 0.2Z & and Salpeter IMF, and 
QSOs (particularly relevant at z ^ 6); Hi) the chemical feedback controlling the prolonged 
transition from Pop III to Pop II stars in the merger-tree model by Schneider; iv) assump- 
tions on the escape fractions of ionizing photons, considered to be independent of the galaxy 
mass and redshift, but scaled to the amount of produced ionizing photons. It then accounts 
for radiative feedback inhibiting star formation in low-mass galaxies. This semi- analytical 
model is determined by only four free parameters: the star formation efficiencies of Pop II 
and Pop III stars, a parameter, r/ esc , related to the escape fraction of ionizing photons emitted 
by Pop II and Pop III stars, and the normalization of the photon mean free path, Ao, set to 
reproduce low-redshift observations of Lyman-limit systems. 

4 Inhomogeneous reionization also produces CMB secondary anisotropies that dominate over the primary 
CMB component for I 4000 and can be detected by upcoming experiments, like the Atacama Cosmology 

Telescope (ACT) or Atacama Large Millimeter/submillimeter Array (ALMA) [28, 29]. 



- 4 - 



A variety of feedback mechanisms can suppress star formation in mini-halos, i.e. halos 
with virial temperatures T v i r < 10 4 K, particularly if their clustering is taken into account 
[43]. It is then possible to assume that stars can form in halos down to a virial temperature of 
10 4 K, consistent with the interpretation of WMAP data [44] (but see also [45]). Even galaxies 
with virial temperature T v i r £ 10 4 K can be significantly affected by radiative feedback during 
the reionization process, as the increase in temperature of the cosmic gas can dramatically 
suppress their formation. 

Based on cosmological simulations of reionization, [46] developed an accurate charac- 
terization of the radiative feedback on low-mass galaxies. This study shows that the effect of 
photoionization is controlled by a single mass scale in both the linear and non-linear regimes. 
The gas fraction within dark matter halos at any given moment is fully specified by the 
current filtering mass, which directly corresponds to the length scale over which baryonic 
perturbations are smoothed in linear theory. The results of this study provide a quantitative 
description of radiative feedback, independently of whether this is physically associated to 
photoevaporative flows or due to accretion suppression. 

Two specific alternative prescriptions for the radiative feedback by these halos have been 
considered: 

i) suppression model - in photoionized regions halos can form stars only if their circular 
velocity exceeds the critical value v cr it = [2A;_BT//um p ] 1 / 2 ; here fi is the mean molecular 
weight, m p is the proton mass, and T is the average temperature of ionized regions, computed 
self-consistently from the multiphase IGM model; 

it) filtering model — the average baryonic mass Mb within halos in photoionized regions 
is a fraction of the universal value = f^/il m , given by the fitting formula M^/M = 
/b/[l + (2 1 / 3 — 1)Mq/M] 3 , where M is the total halo mass and Mq is the total mass of halos 
that on average retain 50% of their gas mass. A good approximation for Mq, is given by 
the linear-theory filtering mass, Mp 3 = (3/a) J Q a da' M 2 J 3 (a') 1 — (a'/a) 1 ^ 2 , where a is the 

cosmic scale factor, Mj = (47r/3)p (irc^/Gp) 3 ^ 2 is the Jeans mass, p is the average total mass 
density of the Universe, and c s is the gas sound speed. 

The model free parameters are constrained by a wide range of observational data. 
Schneider and collaborators reported the best fit choice of the above four parameters for 
these two models that well accomplish a wide set of astronomical observations, such as the 
redshift evolution of Lyman- limit absorption systems, the Gunn- Peterson and electron scat- 
tering optical depths, the cosmic star formation history, and number counts of high redshift 
sources in the NICMOS Hubble Ultra Deep Field. 

The two feedback prescriptions have a noticeable impact on the overall reionization 
history and the relative contribution of different ionizing sources. In fact, although the two 
models predict similar global star formation histories dominated by Pop II stars, the Pop 
III star formation rates have markedly different redshift evolution. Chemical feedback forces 
Pop III stars to live preferentially in the smallest, quasi-unpolluted halos (virial temperature 
^ 10 K), which are those most affected by radiative feedback. In the suppression model, 
where star formation is totally suppressed below i> cr it, Pop III stars disappear at z ~ 6; 
conversely, in the filtering model, where halos suffer a gradual reduction of the available gas 
mass, Pop III stars continue to form at z ;$ 6, though with a declining rate. Since the star 
formation and photoionization rate at these redshifts are observationally well constrained, 
the star formation efficiency and escape fraction of Pop III stars need to be lower in the 
filtering model in order to match the data. Therefore reionization starts at z ,$ 15 in the 
filtering model and only 16% of the volume is reionized at z = 10 (while reionization starts 
at z ~ 20 in the suppression model and it is 85% complete by z = 10). For 6 < z < 7, QSOs, 
Pop II and Pop III give a comparable contribution to the total photo-ionization rate in the 
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filtering model, whereas in the suppression model reionization at z < 7 is driven primarily 
by QSOs, with a smaller contribution from Pop II stars only 

The predicted free electron fraction and gas temperature evolution in the redshift range 
7 < z < 20 is very different for the two feedback models. In particular, in the filtering model 
the gas kinetic temperature is heated above the CMB value only at z ;$ 15. 

The Thomson optical depth, r = J Xe"n e ^Tcdt, can be directly computed for the assumed 
ACDM cosmological model parameters given the ionization histories: tqfo& — 0.1017 and 
t goo — 0.0631 for the suppression and the filtering model, respectively. Note that these values 
are consistent with the Thomson optical depth range derived from WMAP 3yr (7yr) data but 
with ~ la (~ 2a) difference among the two models, leaving a chance of accurately probing 
them with forthcoming CMB anisotropy experiments. 

3.1 Fitting astrophysical histories 

Since filtering and suppression models gives the redshift evolution of ionization fraction 
and electron temperature in a tabular way, as first step we derived these quantities with 
appropriate analytical functions, by means of a specific fitting tool, Igor Pro (v. 6.21) [47]. 
In order to have an accurate parametrization we divided the redshift in bins such to minimize 
the x 2 test given by the fit itself (see Appendix A for details). 

The results found using the complete set of functions specific to each epoch of interest are 
plotted in Fig. 1 for the ionization fraction and in Fig. 2 for the electron temperature. Every 
graph shows the fitting functions (dashed red line) and the corresponding reionization model 
(solid black line). 
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Figure 1. Ionization fraction: comparison between the tabulated data (solid black line) and the fit 
(dashed red line) for the suppression (left panel) and filtering model (right panel). 

The accuracy of the fit can be analyzed by means of the percentile difference among 
theoretical data and fitting functions, namely the ratios between the derived functions and 
the models data, as shown in Fig. 3. From these plots one can see that the difference is 
always < 1%, a bit greater in the filtering model but, in any case, the precision of the fit is 
always very high. 



4 Phenomenological reionization models 

We investigated two phenomenological double peaked reionization models introduced 
by Naselsky &: Chiang in which the Universe was reionized twice at different epochs, the 
first one, late reionization model, at z ~ 10 by Pop III stars and then at z ~ 6 by stars in 
large galaxies, and the second one, early reionization model, at very high redshift (z £ 100), 
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Figure 2. Electron temperature: comparison between the tabulated data (solid black line) and the 
fit (dashed red line) for the suppression (left panel) and filtering model (right panel). 
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Figure 3. Ionization fraction: ratios between the fit functions and the tabulated data for the ionization 
fraction (left panel) and electron temperature (right panel) as function of redshift for the two different 
reionization histories. 

induced for example by decay of unstable particles, followed by reionization mechanisms at 
z ~ 6 as predicted in the standard picture. These models are based on the assumption that 
extra ionizing sources can be described by the efficiency of the ionizing photons production 
coefficient £j which characterizes the ionizing photon production rate: 



^ = £i(z)n b (z)H(z), 



(4.1) 



where is the baryon density and H(z) is the Hubble parameter. 
The late reionization history has been modeled as: 



Ei(z) = e exp 



(z Z re ) 



Az 2 



+ e 1 (l + z)- m Q(z re -z). 



(4.2) 



The two terms in the right hand of the equation describe the two ionizing epochs, with sq, E\ 
and m free parameters of the model, Az the width of the first reionization epoch, the step 
function. The first peak corresponds to a reionization fraction that decreases considerably 
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at z > z re , while the second term characterizes a reionization fraction which is monotonic 
increasing with increasing time (or decreasing redshift). 

For the early reionization history, the authors adopted a representations in which the 
efficiency coefficient has a Gaussian parametric form: 



(4.3) 

again, £, z re and Az are free parameters (clearly different from those of the previous scenario). 

Early Reionization History 
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Figure 4. Dependence of the free parameter £ on the reionization redshift of the early history in 
order to have the adopted value of r. 

Varying the free parameters of the considered model allows us to identify reionization 
histories producing a Thomson optical depth with the desired value. With this approach we 
selected for the late history a set of parameters able to approximately reproduce the optical 
depth found for the suppression model: 



'e = 1.3 • 10 3 

z re = 10 

< P = fa = 5.3 • lCT 6 => t l ~ roFoe = 0.1017. (4.4) 
m = 11.95 
= 0.025z re 



We exploited also the combination of an early reionization history, unable to contribute 
to the actual optical depth having a substantial high reionization redshift, with the filtering 
model such as the global optical depth is exactly the same that figures out from the suppres- 
sion model: 



£i(z) = £exp 



\4 



Az 2 
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'£ = 2.315 

<Z re = 500 => T% t = T E + TGOO = TCF06 = 0.1017. (4.5) 

= 0.025z re 

In Fig. 4 is shown the dependence of the variable £ on the reionization redshift, z re , in 
order to satisfy the above condition (for the considered choice of Az). 

In both astrophysical and phenomenological models we assumed the same cosmological 
parameters defined in the previous section. 

The ionization fraction can be evaluated from the balance between the recombination 
and ionization processes: 

^ = -a rec {T)n b x 2 e + e l (z)(l - x e )H{z) (4.6) 

where 

&rec ~ 4 • 10 ^ 3 (XyiO^-K") ®'®s ^ cm 3 is the recombination coefficient. 

Assuming a curvature term Qk = 0, and so Q\ = 1 — Q, m the Hubble parameter is approxi- 
mated by: 

H(z) = H^£l m (l + zf + n A . (4.7) 

In Fig. 5 we display the time evolution of the ionization fraction for all the models 
considered in this work. In particular, the coupled early and filtering model is plotted for 
three different cases of the reionization redshift. Note that, in order to have a constant 
optical depth, we varied the parameter £ assuming values £ = (2.315,1.031,1.309) when 
z re = (500, 350, 200) respectively, giving a decreasing in the peak of the high redshift region 
of x e , not linear with the reionization redshift. The reason can be found analyzing the redshift 
dependence of the recombination and ionization coefficients that are, respectively: 

Crec 

oc n b oc z , (4.8) 

and 

C^txe^l-ie^oc^ 2 , (4.9) 
so that the higher the lower is x e (see Table 1 for details). 





z re 


Crec 


Cjon 


2.315 
1.031 
1.309 


500 
350 
200 


1.25 • 10 s 
4.29 • 10 7 
8.00 • 10 6 


2.59 • 10 4 
6.75 • 10 3 
3.70 • 10 3 



Table 1. Redshift dependence of recombination and ionization coefficients for decreasing reionization 
redshift and adjusted £ model parameter when joining the early reionization history with the filtering 
model. 

In these prescriptions the ionization fraction and electron temperature are provided 
in an analytical form and connected through the dependence on the temperature of the 
recombination coefficient a rec , that enters in the time evolution equation of the ionization 
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Figure 5. Rcionization fraction: comparison between models. 



fraction. Since the relatively weak dependence, a rec oc T -0,6 , the details of the assumptions 
on the matter temperature are not particularly critical, although, in principle, for the active 
phases out of equilibrium, they could play a non negligible role. For the active phase, we 
assume here, respectively, a temperature profile mimicking the model by Cen in the case of 
late processes and a temperature profile given by the same Gaussian parametric form as in 
Eq. (4.3) but with a peak temperature, T p , as free parameter instead of £ in the case of early 
processes, as adopted by Burigana and collaborators in 2004. We use here T p = 6T0 4 K. When 
the ionizing photons production is negligible and no longer affects the ionization history we 
assume the minimal ionization fraction usually derived in the absence of source terms, thus 
assuring the continuity with the quiescent phases in the evolution of the plasma properties. 

4.1 Comoving fractions of injected photon number and energy density 

For many mechanisms of cosmological reionization, the underlying physical process is 
usually characterized in terms of an additional source of ionizing photons injected in the 
plasma. We link here the parameters of the considered phenomenological histories to the 
comoving fractions of injected photon number density and energy density. Defining the usual 
CMB photon number density uq ~ 399(To/2.7K) 3 , the comoving fractions of photon number 
density injected in the redshift range Z{ > z > Zf is given by: 



An 

n 



dz. 



(4.10) 



The second term, ei(l + z 
easily integrated: 

An 



l Q(z r 



H £j{z)n b {z) 

appearing in the late reionization history can be 



n 



Ein b0 



-dz 



ein b0 



[1 



:i 



(4.11) 



/q no(l + z) l+m nQin 
where n^o = 1.12 • 10~ 5 ^l b h 2 is the current baryon number density. 

For the Gaussian term defining the early model and appearing as first term in the late model, 



- 10 - 



we find the following suitable approximation: 



const • ribo 



An _ 

J Zf n (l + z) 



cxp 



\2l 



Az' 



dz 



const • As. 



no 



1 



7T. 



(4.12) 



where "const" is the early parameter £ or the late parameter sq, respectively, and the as- 
sumption Zf — > — oo and Zi — > oo for the integration limits is made, a good approximation to 
this aim for a peaked Gaussian shape. 

The comoving fraction of injected photon energy density depends on the energy distri- 
bution function of ionizing photons. Given the comoving fraction of injected photon number 
density and assuming the mean ionizing photon energy, E l ° n , here assumed ~ 10 eV for 
numerical estimates, one can compute of the comoving fraction of injected photon energy 
density, Ae/e. In the case of the early history and for the Gaussian term of the late history 
we have: 



Ae 



~ 1.59 • 10" 



AZ7' 



7 



n lOeV \ 2.7 



T, 



:i+z) 



(4.13) 



where for numerical estimates we can assume (1 + 



(1 + 0. 



re ) early /late ■ 



For the second 



term appearing in the late reionization history, after a simple integration we have: 



Ae 



~ 1.60 • 10 



rpiori 
-11 h l 



Ein b0 



lOeV e (m + 1 



■[l-a + Zre) 



-(m+l)l 



(4.14) 



aTi is the current CMB photon energy density. 



where £o — "^o 

In general, we used the D01AJF routine of the NAG Libraries for an accurate numerical 
cross-check of previous analytical formulas and approximations. For the adopted parameters 
we found an agreement better than ~ 0.03%. 



5 Code implementation 

We modified CAMB [48], the cosmological Boltzmann code (see also [49]) for computing 
the angular power spectrum of the anisotropies of the CMB, in order to introduce the ion- 
ization fractions evaluated according to the astrophysical and phenomenological reionization 
models described in previous sections, alternative to the reionization treatment originally 
implemented in CAMB. Of course, the methods described here can be used as guidelines 
for the implementation of any other astrophysical reionization model (see Appendix B for 
further details). 

As a significant step forward with respect to previous analyses, the emphasis of this 
work is posed to the extension to a first detailed characterization of the polarization B-mode 
APS. 

By implementing the source file reionization. f '90, that defines the Reionization module, 
we are able to parametrize the desired reionization history and to supply the corresponding 
ionization fraction as function of redshift [50] . 

Particular care must be taken to the normalization of the quantities occurring in the his- 
tory definition, such as the ionization fraction. In CAMB the reionization fraction is referred 
to the hydrogen, so, when allowing for helium reionization, the global ionization fraction in 
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the case of complete ionization can be greater than one, following the relation: 



x full = 1+ V^L. (5.1) 
n H 

Instead, in the two astrophysical reionization models, the global ionization fraction cor- 
responding to the case of complete ionization is -^f ull ' CF0G = 1.12721 for the suppression 
model, and ^-f ull > G0 ° — 1.12480 for the filtering model. 

To evaluate the total fraction in the case of the phenomenological histories we have 
normalized it with CAMB implementation, such as: 



x C e AMB = ^ = x^ El -^, (5.2) 

where fu is the fraction of baryonic mass in hydrogen. Assuming hydrogen abundance 
fu = 0.76, the global ionization fraction in the case of complete ionization for both histories 
is fh,E = 1.157895. 

Furthermore, it is possible to fix in CAMB the Thomson optical depth parameter, and let 
the code estimates the reionization redshift, or to choose the desired reionization redshift, and 
obtaining from dedicated function the r evaluation. While for the considered astrophysical 
models r is known, we implemented a modified version of this function to derive r according 
to the considered phenomenological model. 

5.1 Astrophysical models and standard CAMB 

In the framework of our adapted version of CAMB code it can be useful to analyze an 
important differences between the various models we described. The crucial characteristic 
resides in the interplay between different physical processes at various epochs which contribute 
to reionize the cosmological plasma. Every mixture of these processes can lead to a distinctive 
scenario, as displayed in Fig. 5. 

Since the standard CAMB traces a reionization that follows the "single peak" evolution 
scheme, it can be interesting to compare it with the suppression and filtering models, inves- 
tigating on their differences in the APS and, in particular, in the sequence of acoustic peaks, 
at various multipoles (see Figs. 6 and 7). 

To this aim we derived the temperature (TT), polarization (EE and BB), and cross- 
correlation (TE) APS of CMB anisotropies for the two astrophysical reionization histories, 
Suppression Reionization Model (CF06) and Filtering Reionization Model (GOO), generically 
denoted as "models" in titles and the legends, and for the original version of the code, denoted 
as "CAMB CF" or "CAMB G", assuming, respectively, an optical depth corresponding to 
the value given by the theoretical model to which we are comparing to. Note also that in the 
right panel of Fig. 7 we display the module of the cross-correlation APS. 

For simplicity, we neglect lensing in this case, and the total Cg plotted represents the 
sum of scalar and tensor contributions. The tensor to scalar ratio of primordial perturbation, 
r, is assumed here equal to 0.1. 

It is also very useful to analyse the relative difference between results obtained with the 
astrophysical models and the original CAMB, as shown in Figs. 8 e 9, defined by the relation: 



riModel r-iCAMB 

1 (5.3) 



l/2(C^ odd + Cf AMB )' 
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Figure 6. Comparison between the models and CAMB in Temperature (left panel) and E-mode 
polarization (right panel) APS. 
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Figure 7. Comparison between the models and CAMB in B-mode polarization (left panel) and 
temperature-polarization cross-correlation (right panel) APS. 



As before, note that in Fig. 9 we display the module of the difference of the cross- 
correlation APS. 

The relative differences found in the case of the suppression model are larger than in 
the case of the filtering prescription, mainly in the polarization and cross-correlation patterns 
of the APS. In addition, the difference is remarkable at low multipoles, in particular at I < 
few tens, i.e. at intermediate and large angular scales, as intuitively expected since we are 
considering relatively late reionization processes. 

6 Sensitivity of CMB measurements and foreground emission 

CMB experimental data are affected by uncertainties due to instrumental noise (crucial 
at high multipoles, £, i.e. small angular scales), cosmic and sampling variance (crucial at 
low £, i.e. large angular scales) and from systematic effects. Also, they are contaminated 
by a significant level of foreground emission of both Galactic and extragalactic origin. In 
polarization, the most critical Galactic foregrounds are certainly synchrotron and thermal 
dust emission, while free-free emission is negligible in polarization and other components, 
like spinning dust and haze, are still poorly known, particularly in polarization. Synchrotron 
emission is the dominant Galactic foreground signal at low frequencies, up to ~ 60GHz, 
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Figure 8. TT (left panel) and EE (right panel) APS: relative differences between the models and 
CAMB. 
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Figure 9. BB (left panel) and TE (right panel) APS: relative differences between the models and 
CAMB. 



where dust emission starts to dominate. External galaxies are critical only at high £, and 
radiogalaxies are likely the most crucial in polarization up frequencies ~ 200GHz, most 
suitable for CMB anisotropy experiments. 

In this section we provide simple recipes aimed at evaluating the levels of sensitivity of 
on-going and future CMB space experiments and the kind of contamination expected from 
foregrounds, and to make us able discuss in the next section on the possibility to distinguish 
between different reionization scenarios in the framework of current and future experiments, 
focussing in particular on the pure polarization E and B-modes. 

6.1 Sensitivity measurements 

The uncertainty on the angular power spectrum is given by the combination of three 
components, Cosmic Variance (CV), Sampling Variance (SV) and Instrumental Noise (N) 
[51]: 



SC e I 2 / Aa 2 
d ~ y f sky (2£ + + NC e W e 

Here f s k y is the sky coverage, A is the surveyed area, a is the instrumental rms noise per 
pixel, N is the total pixel number, Wg is the beam window function that, in the case of a 
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Gaussian symmetric beam, is: 



W e = exp(-£(£+l)al), (6.2) 

being o~b = FWHM/V 8 In 2 the beamwidth which defines the experiment angular resolution. 

For f s ky = 1 the first term in parenthesis defines the cosmic variance, an intrinsic limit 
on the accuracy at which the angular power spectrum of a certain cosmological model defined 
by a suitable set of parameters can be measured with the CMB. It typically dominates the 
uncertainty on the APS at low t because of the small, 11 + 1, number of modes m for each 
t. The second term in parenthesis characterizes the instrumental noise, that never vanishes 
in the case of real experiments. Note also the coupling between experiment sensitivity and 
resolution, the former defining the low £ experimental uncertainty, namely for close to 
unit, the latter determining the exponential loss in sensitivity at angular scales comparable 
with the beamwidth. 

In order to provide concrete estimates of these quantities, we consider Planck LFI and 
HFI channels at v = (70, 100, 143, 217) GHz, and COrE channels at v = (75, 105, 135, 165, 
195, 225) GHz, i.e. at the frequencies particularly suitable for CMB analyses because of the 
combination of good experimental sensitivity and resolution, and of relatively low foreground 
contamination. We adopt here the sensitivities and resolutions summarized in the COrE 
white paper 5 . 



v (GHz) 


FWHM (arcmin) 


ot (/uK • arcmin) 


QPol • arcmin) 


70 


13 


211.2 


298.7 


100 


9.9 


31.3 


44.2 


143 


7.2 


20.1 


33.3 


217 


4.9 


28.5 


49.4 



Table 2. Instrumental sensitivity of Planck experiment. 



v (GHz) 


FWHM (arcmin) 


ax (/uK • arcmin) 


a p i (/uK • arcmin) 


75 


14 


2.73 


4.72 


105 


10 


2.68 


4.63 


135 


7.8 


2.63 


4.55 


165 


6.4 


2.67 


4.61 


195 


5.4 


2.63 


4.54 


225 


4.7 


2.64 


4.57 



Table 3. Instrumental sensitivity of COrE experiment. 



For each of the two projects we computed an overall sensitivity value, weighted over the 
channels, defined by 

^— (6-3) 



a j,tot i °\i 



where j = T, Pol and i states for the sensitivity of each frequency channel, listed in Tables 2 
and 3. FWHM values of 13' and 14' are used to define the overall resolution respectively of 



5 The nominal sensitivity of Planck is slightly better than that adopted here thanks the slightly longer 
extension of the mission with both instruments, and the additional extension with the only LFI. Of course, 
the real sensitivity of the whole mission will have to include also the potential residuals of systematic effects. 
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Planck and COrE in the computation of the beam window function . 

Finally, to improve the signal to noise ratio in the APS sensitivity, especially at high 
multipoles, we will apply a multipole binning of 5% in temperature APS, 15% in TE cross- 
correlation and 30% in polarization APS, both in E and B-modes. 

6.2 Parametrization of residual polarized foreground contamination 

The parametrization of the APS of Galactic thermal dust and synchrotron emission 
adopted in this work is taken from the results of WMAP 3-yrs [52] under the assumption 
that these contributions are uncorrelated 7 , and is expressed by: 

i l±A C { ore = (B.(i//65) 2 ^ + B d (u/65) 2 ^)£ m , (6.4) 

where s and d stands for synchrotron and dust, and the frequency v is expressed in GHz. The 
coefficients characterizing the E and B-modes polarization APS are slightly different, and are 
listed in Table 4. 



E-mode 


b oxy 


P 


m 


B-mode 


B (jjK)* 


P 


m 


sync 


0.36 


-3.0 


-0.6 


sync 


0.3 


-2.8 


-0.6 


dust 


1.00 


1.5 


-0.6 


dust 


0.5 


1.5 


-0.6 



Table 4. Parametrization of E and B mode polarization power spectra of Galactic synchrotron and 
thermal dust emission. 

In the next sections we will parametrize a potential residual from non perfect cleaning 
of CMB maps from Galactic foregrounds simply assuming that a certain fraction of the fore- 
ground signal at map level (or, equivalently, its square at power spectrum level) contaminates 
CMB maps. Of course, one can easily rescale the following results to any fraction of residual 
foreground contamination. The frequency of 70 GHz, i.e. the Planck channel where Galactic 
foreground is expected to be minimum at least at angular scales above ~ one degree, will be 
adopted as reference. 

For what concerns extragalactic source fluctuations [53, 54], we will adopt the recent 
(conservative) estimate of their Poissonian contribution to the APS [55] at 100 GHz 8 assum- 
ing a detection threshold of ~ 0.1 Jy, together with a potential residual coming from an 
uncertainty in the subtraction of this contribution computed assuming a relative uncertainty 
of ~ 10% in the knowledge of their degree of polarization and in the determination of the 
source detection threshold, implying a reduction to ~ 30% of the original level [56]. Except 
at very high multipoles, their residual is likely significantly below that coming from Galactic 
foregrounds. 

7 Results 

In this section we present the angular power spectra resulting from different prescriptions 
of the reionization process and compare them with the sensitivities of on-going and future 

6 In fact, it is possible to smooth maps acquired at higher frequencies with smaller beamwidths to the 
resolution corresponding to the lowest frequency of each experiment. 

7 A more sophisticated treatment of the foreground power spectra takes into account the correlation among 
the various foreground components. 

8 We adopt here a frequency slightly larger than that considered for Galactic foregrounds because at small 
angular scales, where point sources are more critical, the minimum of foreground contamination is likely shifted 
at higher frequencies. 
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space experiments including instrumental noise and fundamental statistical uncertainty. 

Figs. 10, 11, 12, 13 show the temperature, polarization and cross-correlation APS 
up to multipoles t = 2000 for a range of values of the so-called tensor-to-scalar ratio, the 
r parameter (or T/S), which parametrizes the ratio between the amplitudes of primordial 
tensor and scalar perturbations. We also show the estimated cosmic and sampling variance, 
and instrumental noise contributions for the Planck (dotted lines) and COrE (dash-dotted 
lines) experiments, assuming a sky coverage of 80%. 

As expected, the temperature APS does not exhibit a remarkable dependence on the 
T/S parameter, as evident from the comparison between the panels in Fig. 10. Comparing 
the considered reionization models, we observe a difference between their relative power at 
high and low multipoles: in particular, the early plus filtering curves show more power at low 
£ and less at high i. Note that this double peaked reionization history has a high redshift 
phase characterized by a remarkable ionization fraction. 

The E-mode power spectrum is shown in Fig. 11. Again, it is only weakly dependent on 
the T/S parameter, but for early plus filtering model that is widely different from the other 
models up to the second acoustic peak. Note also that, even if suppression and late (double 
peaked) models have the same optical depth and cosmological parameters, their power spectra 
are significantly different in both shape and power. For the late model, the reionization bump 
is slightly shifted to higher multipoles with respect to the other histories. 

The B-mode power spectrum (see Fig. 12), plotted here including also the lensing contri- 
bution, shows the expected linear dependence on T/S where the primordial signal dominates 
over the lensing contribution, which, on the contrary, determines the power at high £, almost 
independently of r. The tensor-to-scalar ratio has also a strong impact on the shape of the 
observable acoustic peaks, in particular the first one is flattened for decreasing values of r, 
because of the relative weight of primordial and lensing signal. Typically, the reionization 
bump is stronger for the suppression and the late (double peaked) models, weaker for the 
filtering and the early plus filtering histories. 

The ideal sensitivity of Planck is enough to detect the primordial B-mode for tensor- 
to-scalar ratios above few x 0.01, in particularly thanks to the information contained in 
the reionization bump and up to the first acoustic peak. The improvement foreseen for 
an experiment with a sensitivity like COrE could allow to reveal the primordial B-mode 
polarization signal down to r ~ 0.001 (or even lower). 

The ultimate limitation comes from foregrounds. In the case of the B-mode, we show an 
estimate of the contamination by Galactic synchrotron and thermal dust polarized emission, 
and by extragalactic point source fluctuations, parametrized as described in the previous 
section. In all cases, the extragalactic signal is never dominant except at very high multi- 
poles, but still remaining below the contribution by lensing, while Galactic foreground may 
significantly contaminate the CMB measure, especially at low and intermediate multipole. 

The set of panels in Fig. 13 presents the temperature-polarization cross-correlation APS, 
plotted in absolute value. The plus sign at the top of each panel denote positive correlation. 
Again, there is only a weak dependence on the T/S parameter, and a substantial difference 
between the considered evolutionary models, in particular for the early plus filtering history 
up to t ~ 400. 

7.1 Comparison between models 

In order to understand if different models can be distinguished, we analyze the relative 
differences between two models and compare them with experimental sensitivity and fore- 
ground residual estimates (properly normalized). 
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Figure 10. TT APS for the reionization histories: suppression, filtering, late double peaked, early 
plus filtering for different values of tensor to scalar ratio (see plot legend and text for details). 



Since we are considering here several models, it is useful to define a reference model 
to which divide their differences as well as the experimental sensitivity and foreground lim- 
itation. We assume here a ACDM model with the same cosmological parameters adopted 
in the models under investigation, and an optical depth fixed by the suppression history 
( T Rif = t cf)- We compute its power spectra using the standard CAMB. 

Note that, in principle, each model could be normalized in order to match available 
CMB data in the desired range of multipoles. Currently, temperature data play the major 
role in this respect. The overall normalization of the APS of a given model is related to 
the amplitude of initial perturbations, or, equivalently, to the density contrast at a reference 
scale, typically assumed at 8h _1 Mpc, i.e. the parameter as, which is better determined when 
CMB data are combined with galaxy surveys [57]. In practice, this involves a multiplicative 
factor of the APS. According to this choice, the relative differences between two models could 
be more or less prominent at different multipole ranges. In order to make our comparison 
between models dependent only on the APS shape and not on the normalization adopted 
for each model, before of computing their relative differences, we renormalize each model to 
make its TT APS averaged over I equal to that of the reference model. 
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Figure 11. EE APS for the reionization histories: suppression, filtering, late double peaked, early 
plus filtering for different values of tensor to scalar ratio r (see plot legend and text for details) . 



We report in Figs. 14, 15, 16, 17 the relative differences between models for the 
temperature anisotropies, polarization EE and BB APS, and the temperature-polarization 
cross-correlation, respectively. Each panel shows the comparison between two models, nor- 
malized to the reference model defined above, for a wide set of the T/S parameter. 

As anticipated, the tensor-to-scalar ratio does not affect significantly the temperature 
anisotropies, and for this reason in Fig. 14 the curves appear almost superimposed. 
In addition, at £ £ 100 — 300 the differences tend to be approximately null, except for the 
comparison of the early plus filtering model with all the others, because of their dramatic 
difference at high redshifts. 

Considering the whole multipole range, the largest differences appear again in the comparison 
of the early plus filtering model with all the other (later) histories, achieving a maximum level 
of ~ 70% at £ ~ 10 - 15. 

Note that we can discriminate only early processes from other (later) models using only 
the CMB TT APS. This can not be significantly improved with the future generation of 
experiments, the limitation coming essentially from cosmic variance. Thus, polarization data 
are crucial. 



- 19 - 



Total BB & Lensing, r=0.1, fsky=0.8, bin 30% 




Total BB & Lensing, r=0.3, fsky=0.8, bin 30% 



Suppression, t = 0.1017 
Planck CV+N 
COrE CV+N 
Filtering, t = 0.0631 
Planck CV+N 
COrE CV+N 

Synchrotron, z/ cmb = 70 GHz 
Radiosources 




Total BB & Lensing, r=0.003, fsky = 0.8, bin 30% 

• • ' 77 




Late Double Peaked, t = 0.1017 
Planck CV+N 
COrE CV+N 

Early & Eiltering, t = 0.1017 

Planck CV+N 

COrE CV+N 

Dust, iv cmb = 70 GHz 

Radiosources 30% 



Figure 12. BB APS for the reionization histories: suppression, filtering, late double peaked, early 
plus filtering for different values of tensor to scalar ratio r. Galactic synchrotron (dash-dotted blue 
line) and thermal dust (dash-triplc-dotted cyan line) polarized emission, and extragalactic point source 
fluctuations (solid orange line) and their potential residual (solid violet line) as described in previous 
section (see plot legend and text for details). 



The comparison between the E-mode polarization power spectra is more interesting (see 
Fig. 15). As evident, in particular in the three last panels in Fig. 15 where the early plus 
filtering scenario is compared with the others, the differences related to the r parameter 
emerge clearly. The greater is r the greater are the relative differences. The large difference 
at £ ~ 200 is due to the first early reionization phase, absent in the other models. 

In general, remarkable relative differences appear between all the models at t ~ 5 — 25 
or even up to t ~ 100 when comparing the early plus filtering scenario with the others. 
They are connected to the details of the ionization history at lower redshifts. We also plot a 
potential residual contamination from Galactic foregrounds. Note that a foreground removal 
at a few per cent level of accuracy at map level makes astrophysical contamination below the 
sensitivity limitation of on-going and future space experiments. The difference in nominal 
capability of Planck and COrE is not so remarkable in this respect, although, in practice, a 
significant improvement in sensitivity and frequency coverage clearly will make much more 
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Figure 13. TE APS for the reionization histories: suppression, filtering, late double peaked, early 
plus filtering (see plot legend and text for details). Here is plotted the absolute value of the TE cross 
correlation, the "blue plus" indicates where the Ce are positive and each cusp is an inversion point in 
the sign of the Ci themselves. 



robust and accurate the foreground subtraction process. 

More complex is the case of the B-mode polarization power spectra, reported in Fig. 16, 
where, for completeness, we display also an estimate of the potential residual from extragalac- 
tic point source fluctuations which is always negligible in comparison with the other sources 
of limitation. As before, the early plus filtering model largely differ from the others. Note 
that for the B-mode the relative differences between the various models remarkably depend 
on r in all cases. In this representation, the Planck sensitivity strikingly depends on r at low 
and intermediate £, because of the instrumental noise limitation, while that of COrE is almost 
independent on r, being essentially cosmic variance limited where the lensing contribution 
does not dominate. The sensitivity of COrE instruments is such that we could be able to 
distinguish between suppression and filtering (or late) models in a certain range of multipoles 
(around I ~ 10) even for r slightly larger than ~ 10~ 3 , while the early plus filtering history 
can be distinguished from all the other histories on a wider multipole region (or for even 
lower values of r). The impact of residuals of Galactic polarized foregrounds clearly increases 
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for decreasing tensor-to-scalar ratio, as expected in the comparison of relative differences 
between models. Only for r ^ 0.1 a foreground subtraction at 3% accuracy level in the map 
is enough to discriminate between each couple of the considered models, while the early plus 
filtering history can be distinguished from all the other histories even with a less accurate 
foreground subtraction, thanks to its prominent differences at intermediate multipoles. 

Finally, relative differences in the temperature-polarization cross-correlation APS are 
also very large on very wide ranges of multipole when comparing the early plus filtering 
history with all the other models (see Fig. 17), filtering and late models show remarkable 
differences (larger than 80%) at low multipoles, while the suppression model differ from the 
late and filtering models only on a very small range of multipoles around £ ~ 10. The 
"spikes" appearing at high multipoles are due to little shifts of the multipoles corresponding 
to the change of sign of the cross-correlation spectra for the different considered models. The 
difference in nominal capability of Planck and COrE is not so remarkable for the TE mode. 

8 Conclusion 

The inclusion of astrophysically motivated ionization and thermal histories in numerical 
codes is crucial for the accurate prediction of the features induced in the CMB, for constrain- 
ing reionization models with CMB data, and to exploit current and future high quality CMB 
data with great versatility to accurately extract cosmological information. 

We have implemented a modified version of CAMB, the Cosmological Boltzmann code 
for computing the APS of the anisotropies of the CMB, to introduce the hydrogen and helium 
ionization fractions predicted in two astrophysical reionization models, i.e. suppression and 
filtering model, in two classes of phenomenological reionization histories, involving late or 
early reionization, and in their combination, as alternative to the original implementation of 
reionization in the CAMB code, and beyond the simple r-parametrization. For astrophysical 
models, we provide also suitable analytical descriptions of the ionization and thermal his- 
tories that can be ingested in any numerical code aimed at computing CMB features. We 
compared the results obtained for these models for all the non-vanishing (in the assumed 
scenarios) modes of the CMB APS. As a significant step forward with respect to previous 
analyses, the emphasis has been posed here to the extension to a first detailed characteri- 
zation of the polarization B-mode APS. Its amplitude and shape depends, in particular, on 
the tensor-to-scalar ratio, r, and on the reionization history, thus an accurate modeling of 
the reionization process will have implications for the precise determination of r or to set 
more precise constraints on it through the joint analysis of E and B-mode polarization data 
available in the next future and from a mission of next generation. This is particularly crucial 
in the case of low values of r, because of the contribution to B-mode coming from lensing 
that competes with the primordial B-mode at intermediate multipoles. Taking into account 
also the limitation from potential residuals of astrophysical foregrounds, we discussed the 
capability of next data to disentangle between different reionization scenarios in a wide range 
of tensor-to-scalar ratios. 
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Figure 14. Relative differences between the astrophysical and phcnomenological reionization histories 
in temperature power spectrum for all values of the r parameter assumed in this work (see plot legend 
and text for details). is the adopted normalization power spectrum. 

Appendices 

A Fitting ionization and thermal histories 

The considered astrophysical reionization models, as well as others published in the liter- 
ature, provide ionization and thermal histories in tabulated form. The cosmological analysis 
of CMB anisotropy data largely relies on Boltzmann codes for computing the angular power 
spectrum in temperature, polarization, and cross-correlation modes under general conditions. 
The inclusion in such codes of reionization histories beyond the simplistic phenomenological 
approximations already implemented in the publicly available codes, allows to achieve more 
accurate predictions, of particular interest for the analysis of polarization data. Having func- 
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Figure 15. Relative differences between the astrophysical and phcnomcnological reionization histories 
in polarization EE-mode power spectrum for all values of the r parameter assumed in this work. C e 
is the adopted normalization power spectrum. Potential residuals of Galactic foregrounds are also 
shown (see plot legend and text for details). 



tional descriptions of the evolution of the ionization fraction allows to speed-up computation 
and improve code versatility with respect to the use of interpolation of tabulated grids. Al- 
though in this context only the ionization history is relevant, we report here for completeness 
also the results concerning the thermal history. 

The software we used to fit to the suppression and filtering model is Igor Pro (v. 6.21 ), an 
integrated program for visualizing, analyzing, transforming and presenting experimental data, 
such as curve-fitting, Fourier transforms, smoothing, statistics, and other data analysis, image 
display and processing, by a combination of graphical and command-line user interface. We 
report here all the functional forms we found particularly suitable to represent the considered 
ionization and thermal histories, since they could be useful as guidelines to fit other kinds of 
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Figure 16. Relative differences between the astrophysical and phcnomenological reionization histories 
in polarization B-mode power spectrum for all values of the r parameter assumed in this work. 

is the adopted normalization power spectrum. Potential residuals of Galactic foregrounds and 
extragalactic point source fluctuations are also shown (see plot legend and text for details). 

reionization history. More details about the usage of this software to the current aims and 
the complete list of parameters of the below functions are given in [50] . 

A.l Fitting functions for the reionization histories 

In the case of the suppression model the redshift intervals and the relative analytic 
functions are: 

- Interval z = [0, 3.8] - Polinomial Function of 6° order: 

Xre = ao + a\z + a 2 z 2 + a 3 z 3 + a 4 z 4 + a 5 z 5 + a G z 6 . 
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Figure 17. Relative differences between the astrophysical and phcnomenological reionization histories 
in temperature polarization cross correlation power spectrum for all values of the r parameter assumed 
in this work (see plot legend and text for details) . Cf 1 ^ is the adopted normalization power spectrum. 



- Interval z = [3.8, 6] - Polinomial Function of 5° order: 

Xre = b + b x z + b 2 z 2 + b 3 z 3 + 6 4 z 4 + b 5 z 5 . 

- Interval z = [6,9] - Polinomial Function of 5° order: 

Xre = Co + C\Z + C 2 Z 2 + C 3 Z 3 + C 4 Z 4 + C 5 Z 5 . 
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- Interval z = 



[9, 12.4] - Log-Normal Function: 

f (X 

Xre = lo + hew - - 



Interval z = [12.4, 14.2] - Sigmoidal Function: 

di 



Xre = d + 



l + exp(^ 



Interval z = [14.2, 16.8] - Hill Equation: 

Xre = e + (ei - e ) 



yt'2 



z e 2 + e% 



'-■2 



Interval z = [16.8, 18] - Decaying Exponential Function (expXOffset) 

'h s - z s 



Xre = h + /iiexp 



Interval z = [18, 20] - Hill Equation: 



Xre = 10 + (h - *o) 



z^ + i% 



Interval z = [20, 20.2] - Linear Function: 



Xre = m + m\Z. 



Interval z = [20.2,23] - Decaying Exponential Function (expXOffset) 



Xre = fo + /iexp 



h 



Interval z = [23, 30] - Log-Normal Function: 



Xre = 9o + ffiexp 



ln(z/g 2 ) 



93 



When considering the filtering model the fitting functions are: 
- Interval z = [0, 3.8] - Polinomial Function of 9° order: 



Xre = o,q + a\z + Q2Z 2 + a^z 3 + a^z 4 + as^ 5 + agz 6 + a-jz 1 + agz 8 + agz 9 . 



Interval z = [3.8, 6] - Log-Normal Function: 



Xre = io + nexp 



ln(z/i 2 ) 



>3 



Interval z = [6, 6.2] - Linear Function: 



Xre = m + miz. 



Interval z = [6.2, 9] - Hill Equation: 



Xre — h + (h — h) 



Interval z = [9, 11.6] - Sigmoidal Function: 



Xre = b + 



bi 



l + exp(^ 



Interval z = [11.6, 13] - Power Function: 

Xre = C + CiZ C2 . 

Interval z = [13, 15] - Power Function: 

Xre =h + U\Z H% . 



Interval z = [15, 17] - Decaying Exponential Function (expXOffset): 

'd 3 - z" 



Xre = d + diexp 



Interval z = [17, 19.6] - Decaying Exponential Function (expXOffset): 

e 3 - z 



Xre = e + eiexp 



C 3 2 
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Interval z = [19.6,22.2] - Decaying Exponential Function (expXOffset): 

'h-z 



Xre = fo + /iexp 



h 



Interval z = [22.2,30] - Decaying Exponential Function (expXOffset): 

'93 ~ z 



Xre = 9o + 9iexp 



92 



A. 2 Fitting functions for the temperature histories 

In the case of the electron temperature, the fitting functions for the suppression model 

are: 

- Interval z = [0, 3.8] - Polinomial Function of 8° order: 

T re = (iq + a\z + d2Z 2 + fl3-z 3 + a^z 4 + (15Z 5 + clqz® + cijz 7 + agz 8 . 



Interval z = [3.8, 5.8] - Sigmoidal Function: 



T re — bo + 



l + exp(^ 



Interval z = [5.8, 13] - Hill Equation: 



T re = c + (ci -c ) zC2 +c . 2 



Interval z = [13, 16] - Hill Equation: 

T re = h + {h — lo)- 



Interval z = [16, 16.8] - Decaying Exponential Function (expXOffset): 

'd 3 - z" 



T re = do + eqexp 



- Interval z = [16.8, 18.6] - Hill Equation: 



T re = e + (ei - e ) 



z^ + e e 3 2 
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Interval z = [18.6,20.2] - Log-Normal Function: 



T re = ho + /iiexp 



Hz/h 2 ) 



Interval z = [20.2,21] - Sigmoidal Function: 



T re — io + 



11 



1 + exp 



is 



Interval z = [21,22] - Log-Normal Function: 



T re =9o + 9i exp 



ln(z/g 2 ) 



.93 



Interval z = [22, 30] - Log-Normal Function: 



T re = fo + /iexp 



h 



Finally, the fitting functions for the filtering model are: 

- Interval z = [0, 3.8] - Polinomial Function of 8° order: 

T re = + oi-z + a 2 z 2 + (13Z 3 + a^z 4 + az,z 5 + a^z 6 + ajz 7 + agz 8 

- Interval z = [3.8, 6.2] - Sigmoidal Function: 



T re — bo + 



1 + exp ' ^£ 



Interval z = [6.2, 10] - Log-Normal Function: 



T re = C + Ciexp 



V c 3 



Interval z = [10, 11] - Decaying Exponential Function (expXOffset): 

'h 3 - z" 



T re = ho + /iiexp 



Interval z = [11, 15] - Sigmoidal Function: 
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T re — do + 



di 



1 + exp W 



Interval z = [15, 18] - Sigmoidal Function: 

r re = eo + 



ei 



1 + cxp 



Interval z = [18, 21] - Log-Normal Function: 



T re = fo + /iexp 



V /3 



Interval z = [21,30] - Log-Normal Function: 



T re = 9o + 5iexp 



ln(«/ff2) 



9-4 



B Code implementation in CAMB 

We give here some details on the routines we implemented in CAMB to make it able 
to reproduce the considered histories for the ionization fraction and electron temperature, 
summarizing the main innovations included in the modules of interest. 

B.l Subroutine modifications in Reionization module 

The first improvement concern the type ReionizationParams with the inclusion of a new 
string variable, history, through which the user can discriminate between the models, stored 
in the settings parameters file params.ini. 

The main function of the original module, Reionization jxe, has been now written for 
each history, Reionization JxeC F , ReionizationjxeG , Reionization _xeL, Reionization jxeE . 
All of them retrieve the analytical ionization fraction for each redshift bin of interest. 
With this approach is not necessary to parametrize Xe in terms of the variable Window- 
VarMid: 



y = {l + zf'\ (B.l) 

where the exponent is set with the constant Reionization_zexp. 

Some initial parameter values have been revisited too, such as the maximum redshift 
at which \e varies, fixed to 700 instead of 40 to take into account high redshift reionization 
phases, necessary for the early history, and the corresponding initial scale factor astart, 
inversely proportional to the redshift: 



1 



l + z 



(B.2) 



In the function Reionization-timesteps the minimum number of time steps between 
taujstart and tau-complete, the relevant times for the reionization process, has been incre- 
mented to 1000, while for the implementation of the adopted ionization histories the functions 



- 31 - 



listed below are no longer necessary: 



- Reionization jdoptdepth^z) , the subroutine that expresses the integral optical depth in terms 
of the scale factor, 

- Reionization_GetOptDepth(Reion, ReionHist) , the routine which evaluates the integral of 
the optical depth in the redshift interval (0, z™), 

- Reionization_zreFromOptDepth(Reion, ReionHist), a general routine to find the z re pa- 
rameter given optical depth, 

- Reionization_SetFromOptDepth(Reion, ReionHist), the subroutine that calculates the 
redshift of reionization. 

This set of function, in fact, is related to the optical depth definition, i.e. 



as computed in the standard CAMB, while for the histories under examination we have imple- 
mented specific codes to evaluate it, both in our modified CAMB version and as independent 
codes. 

B.2 Subroutine modifications in ThermoData module 

The Thermodata module, implemented in modules. f90 source file, contains the subrou- 
tine inithermo(taumin, taumax), which evaluates the unperturbed baryon temperature and 
ionization fraction as function of time. If there is reionization, the function discriminates be- 
tween the models, smoothly increases Xe to the requested value and sets the actual joptAepth 
to the value imposed by the corresponding model. 
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